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Abstract 

We consider an Unruh-DeWitt detector interacting with a massless 
Dirac field. Assuming that the detector is moving along an hyperbolic 
trajectory, we modeled the effects of fluctuations in the event horizon 
using a Dirac equation with random coefficients. First, we develop the 
perturbation theory for the fermionic field in a random media. Further 
we evaluate corrections due to the randomness in the response function 
associated to different model detectors. 
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1 Introduction 

In a previous paper [I], the spectral density associated to a scalar field 
assuming event horizon fluctuations was presented. In particular, it was 
considered a two-level uniform accelerated system interacting with a massless 
scalar field. Assuming the two-level system prepared in the ground state and 
the field in the Minkowski vacuum state, the main result of the paper was to 
show how the transition rates are modified by event horizon fluctuations. For 
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modeling the event horizon fluctuation it was assumed that the scalar field 
satisfies a massless Klein-Gordon equation with random coefficient. It was 
found a modified response function with the same local temperature found in 
the non-fluctuating case. The correction due to the event horizon fluctuations 
has a Fermi-Dirac factor. A similar result found previously by Takagi [2] in a 
quite different situation. For a careful discussion for the apparent inversion 
of statistics see for example the Ref. [3J. 

The main motivations for such kind of studies follows the ideas present 
in Ref. [3]. In this work it was proposed in condensed matter physics an 
analog model for quantum gravity effects which in turn is based on results 
obtained by Ford and collaborators E] arid Hu and Shiokawa [TJ. Two 
general features of waves propagating in random fluids were considered. 
First, acoustic perturbation in a fluid define discontinuity surfaces that 
provide a causal structure with sound cones. Second, propagation of acoustic 
excitations in random media are generally described by wave equations 
with random speed of sound [H [9], [101 [TTJ. A quantum scalar field theory 
associated to acoustic waves was analyzed in a situation where the velocity of 
propagation of the acoustic wave and consequently the sound cone fluctuates. 
In this approach a stochastic ensemble of fluctuating geometries was assumed. 

In this paper we follow the ideas present in Ref. [TJ. Basically we study 
how a particular model for the fluctuations of a black-hole event horizon 
can affect the transition rate of a two-level system [121 EBJ. Let us consider 
the line element of a four-dimensional Schwarzschild space-time describing a 
non-rotating uncharged black-hole of mass M: 

,h»=(l-^)*»-(l-^)"V-r»*l» (1) 

where dQ 2 is the metric of a unit 2-sphere. Close to the horizon r « 2M, 
the line-element given by Eq. flTJ can be written as 

ds 2 = (jj^J dt 2 - dp 2 - AMdtt 2 , (2) 

where p(r) = ^/8M(r — 2M) and the quantity AMdVt 2 describes the line 
element of a 2-sphere of radius AM. In turn, the other contribution can 
be identified with the line element of the two-dimensional Rindler edge by 
setting t = AMar and p = — , for < p < oo and — oo < t < oo. 
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We can summarize our approach saying that one model to take into 
account the effects of black hole event horizons fluctuations is the following: 
we assume that event horizon fluctuations are modeled basically by a 
centered, stationary and Gaussian random processes. Under these conditions, 
we obtain a random differential equation for the Dirac field that cannot be 
solve exactly Therefore, we implement a perturbation theory, similar to the 
one developed in Ref. [2], in our case associated to a massless Dirac field. 

The organization of the paper is as follows: In section [2] we discuss the 
perturbation theory for the fermionic field in a disordered medium. In section 
[3] we present the modified positive Wightman function due to the fluctuating 
horizon. In section H] we study the modifications to transition probabilities of 
a Unruh-DeWitt detector due the event horizon fluctuations. Finally, section 
[5] contains our conclusions. To simplify the calculations we assume the units 
to be such that K — c — ks — 1. 



2 Dirac field in disordered medium 

We begin by writing a random Dirac equation 

d 

i T °(l + /i(x))— -i^Vi-m 

Here, \I/(£, x) represents a fermionic field, 7°,7* are the Dirac matrices and 
/x(x) is a dimensionless random function of the spatial coordinates. We 
consider a zero-mean random function, i.e., 

<A*(x)>„ = (4) 
and for simplicity we suppose white-noise correlations, i.e., 

(Mx)/i(x')>M = ^(x-x')- (5) 

The symbol denote ensemble averaged of noise realizations and a 

represents the strength of the noise. We also suppose that the noise is 
Gaussian distributed. These last statements are pragmatic one. This is 
one of the simplest models one can choose which can exhibit light cone 
fluctuations. We will also see in next sections that this choice implies in 
a loss of stationarity in the Wightman function to the Rindler noise and the 
Wiener- Khint chine theorem is no longer valid [To'] . 



tf(t,x) = 0. (3) 
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To proceed, let us define the Fourier transforms to the Dirac field and 
also to the noise: 

I" rlk 

Mx) = J j^y^^^k). (7) 

Similar to Ref. [H] we can split Eq. Q in a free and perturbed 
part (Lq + Li)ty = 0, where Lq and L\ are matrices in the coordinates 
or momentum spaces. We have 

L = i7°^-i7 i Vi-m (8) 

and 

Li = t7VW^. (9) 
In the momentum space these operators are written as 

Co = [^u + 7% - m] 5(k - k') (10) 

and 

d = 7 °w/i(k-k'). (11) 

Thus, it is possible to define the full (operator valued) fermionic Green 
function, S = (Cq + Ci) -1 . Assuming a weak noise, a perturbative expansion 
for S is given by 

S — So — SqjCiSq + SqCxSqCxSq — ••• (12) 

where So is just the inverse of the free operator and it is Fourier transformed 
by: 

Sb&x) = y d*ke- ik ^So(u,k). (13) 

In the massless case, So (a;, k) is given by the well known expression: 

a f in 7oW + 7^< 

5 (w,k) = g 14 

ur — k z 
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After performing the average over the fluctuations, which can be seen 
also as an average in the ensemble of geometries in Eq. (1121) . we can write 
down the connected two-point Green function related to the Dirac field as a 
Dyson equation 

(S) li = S -S E(S) li , (15) 

where £ is the irreducible self-energy. The first contribution given by 
S^(t,x) is of second order in random function fi. In the Fourier 
representation it assumes the form given by: 

5 (1) (t,x) = J 0- 4 S^(u,k)e- ik ^'\ (16) 

where 

S"V,k) - (17) 
In the above equation the self-energy E(u;, k) is given by 

E ( M ,k) = /jk- ^:^ . (is) 

Only the first part of the integral in Eq. (fTBl contributes giving 2m 2 uj 2 ^ . 
A detailed calculation given in appendix A shows that the correction to the 
positive frequency Wightman function is given by 

S (1) M = £rfh + ^ (T 1 + 7 2 + 7 3 ) h (19) 

which 

I x = -3684At 6 + 7539At 4 Ax 2 - 8328At 2 Ax 4 - 1005Ax 6 (20) 
+ze(18504At 5 - 24732At 3 Ax 2 + 9354AtAx 4 ) 

and 

I 2 = -1660 At 5 Ax - 2024At 3 Ax 3 - 156AtAx 5 (21) 
+k(6620At 4 Ax + 5056 At 2 Ax 3 + 12Ax 5 ). 
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We have used for Ii and I2 an first order expansion: (x + a) n ~ x 11 + nx n ~ l a 
and D has been defined by D = (At — ie) 2 — Ax 2 . 

We are now able to apply this result to the calculation of the response 
function of the Unruh-DeWitt detector. This will be done in next sections. 

3 Wightman functions and Rindler noise 

Before consider the Unruh-DeWitt detector we need to understand how 
the Wightman function and the Rindler noise changes if we take into 
account event horizons fluctuations. First, it is known that for an uniformly 
accelerated reference frame with proper acceleration a, Rindler's coordinates 
are the appropriate coordinates to be used. The relation to the Minkowskian 
coordinates is given by 

t = a -1 sinh ar 

x = a -1 cosh ar (22) 

where r is the observer's proper time. In the observer's proper reference 
frame it is necessary to Fermi- Walker transport the Dirac field due its non- 
stationarity [2]. We define a proper spinor in the following form 

%) = S t V(x(t)), (23) 

where the corresponding matrix S T takes care of the Fermi- Walker transport. 
Here the observer is, by definition, a Rindler observer, in a given time r 
related to the laboratory frame we define the boosts A T , which the non- 
vanishing components are given by 

(A T )J = (A T )} = coshar, 
(Ar)® = (A T )l = -sinhar, 

(K)l = (A T )l = I- (24) 
So the corresponding matrix of the transformation of the spinor is given by 

S T = e i«*wn = cosh(ar/2) + 7071 sinh(ar/2). (25) 
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The Wightman function will be also Fermi- Walker transported due its 
definition in terms of the Dirac Field. The "proper" Wightman function in 
the accelerated reference frame is called the Rindler noise for the Dirac field 
and it is defined by: 

S(t, t') = S t (S + (x(t), x(t')))^. (26) 

Changing the variables £ = r — r' and 77 = r + r' the total Green function 
is given by 

S(^ V )=S^^r ] )+S^^r ] ). (27) 

The Rindler noise associated to the free case (non-random Dirac 
equations) is given by 



S®{Z,V) = —J-Y \l° cosh (22) +y sinh (ZL)] S ^f) , (28) 
/; 2tt 2 V27 I 1 \2J 1 V 2 7J S inh 4 (f -ze) V ! 

which is the result already known in literature. Using Eqs. (ll9j) . (|20|) and ( l2Tj) 
we find that the Rindler noise, taking into account fluctuations of the event 
horizon is given by: 



sinh — iei) 
sinh(^) 



'sinh 7 ('^r — i^2i 
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where some terms were absorbed by the e« factor. The strength noises 
were rewritten by of (7/) = —a 2 F(rj) and of (7/) = —a 2 H{rj) which are 
positive functions and grow with 77, while the functions F{rj) and H{rj) are 



given by Eqs. flBTT) and (JHSD in Appendix B This result shows how the 



fluctuating event horizon modifies the Rindler noise associated to the Dirac 
field. In the next section we will apply these results and the results found 
in previous section to find the response function associated to the Unruh- 
DeWitt detector. 
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4 Unruh-DeWitt detectors for Dirac 
particles 

In this section we use our results obtained with the Dirac field to study 
the response and transition probabilities associated with the Unruh-DeWitt 
detector [TO] and [TTJ . We assume an uniform accelerated detector coupled 
to a massless Dirac field ^> with the following interacting Hamiltonian 

H int = c im (T)^(T)^(r). (30) 

In the Eq. ( 1301) x M (r) is the world line of the two level system with respect 
to the proper time r, m(r) is the monopole moment operator and \1/ denotes 
a Dirac conjugate to the field and c± is the coupling constant between the 
detector and the Dirac field. 

We consider the detector as a point-like object with a two energy levels 
structure given by u g < u e and eigenstates \g) and |e) respectively [IS] and 
[15] . The gap energy between both states would be E = u e — u g . Defining 
the initial state of the system at r = as |tj) = \g) £g> \ipi) and a final state at 
time r by |tj) = |e) <S> \ ipj) we can apply perturbation theory to compute the 
probability of transition of the Unruh-DeWitt detector. Let us assume that 
the Dirac field is prepared in the Minkowski vacuum state \0,M). In first 
order perturbation theory, we get that the probability of transition P(r, 0) 
is given by 

P(r,0) =c 2 1 \(e\m(0)\g}\ 2 R(E 1 T 1 0), (31) 

where the detector selectivity is given by | (e|?Ti(0) | (?) | 2 . Now we will only 
concentrate in the response function which is given by 

R(E, t, 0) = r dr' [ T dr"e- iE ^'- T "^ (32) 
Jo Jo 

x (0, M|^(r / )^(r')^(r' / )^(r' / )|0, M). 

Then, in this case we have to calculate a four-point correlation function. Due 
to the Gaussian nature of the noise this four-point correlation function can 
be compute in terms of a product of two-point correlation functions. In this 
case, one can find [19] : 
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(0,M|W(r , )*(/)W(r")*(r")|0,M) ^Tr[(S + (r' } r")) 2 ] : (33) 

where Tr is the trace and S' + (r', r") = (0, M|^(r / )^(r")|0, M) is the positive 
frequency Wightman function. Introducing again the variables £ = r' — r" 
and 7/ = t' + r" the response function may be written as [20] 

J2(E, r) = \f dH f T dr/e-^Tr[(5 + (£, r/)) 2 ]. (34) 

After averaging over the fluctuations, we can introduce Eq. ( 1271) into 
Eq. (1341) to get the one-loop correction to the response function. Hence, 
we expand the squared Wightman function in Eq. (1331) resulting in four 
contributions to the response function, where omitting the arguments of the 
Wightman functions it can be expressed by: 

Tr((S + ) 2 ) = Tr(S 2 + S S {1) + S {1) S + (S {1) ) 2 ). (35) 

Substituting the Eq. ([35]) in Eq. flM} we can write R(E, r) = R (E, r) + 
R\(E, r), where 



i r r /• 2r -l?l 
R (E,r) = -J d£j^ d V e-^Tr(S 2 ). (36) 



and 



1 fT f2T-\£\ 

R 1 (E,t) = - dU d V e-^Tr(S S^+S^S ). (37) 

1 J~T J\i\ 



In the asymptotic limit the first term R (E,t) gives the known thermal 
contribution. We show only the result because it can be calculated following 
the same steps done for the correction which we will show later. The main 
aspect is that in the free case (non-random Dirac equation) the probability of 
excitation, i.e., E > (modulo the selectivity) of the Unruh-DeWitt detector 
presents a Bose-Einstein factor [19]. We get lim T ^. 00 Ro(E , r) = W (E,t) 
where 

1 Et 

Wo(E, r) = ^^T-T(4« 4 + 5£V + E% (38) 

Note that the second and third terms in Eq. ( |35l) are the main corrections 
up to a 2 order, see Eqs. ( 1281 and ( 1291) . and the last term is of the order a 4 . 
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So, taking the limit e — > 0, the leading contribution due to horizon events 
fluctuations is given by: 



" ;2(r ?Ui 2(7?) > 




2^\2) / L, ' sinh 9 (¥' 



where we defined <j'i {ff) = <j\ (77) cosh (^) and a'^iv) = a 2(v) s i nn (ip)- 
The integration over 77 it's easier to be done using algebraic computational 
programs. Let us write these integrations in terms of two functions /i(r, £) 
and /2(t, £). Now we can obtain the expression: 

fll(B , T) = 4 r d(e ->* fAwwM) (40) 



To perform these integrals it is convenient to replace \ = at; and express 
each integral as 



1 f aT _E 1 

- / dxe~^ x h(r, x ) = - \ dxe~ l ^ x h(r, x) 

./ — rvr 



1 

a 



CO 



dxe l « x h(r,x)+ / dxe ^ x /i(r,x) 

(41 



where h(r,x) are the terms in the curly brackets from Eq. (j4"0l) . Replacing 
X — > —X m the first integral into squared brackets, we can simplify all 
brackets to a single integral given by: 

^• T )=~(i)'r* sm (f) ft(T ' x) - <42) 

The nature of this term is associated with the switching on and off of 
the interaction between the detector and the background field. It can be 
neglected if we assume large proper time intervals. However, for short proper 
times intervals it must be taken into account. 

We need now to evaluate the infinity range integrals in Eq. (14"T|) . These 
integrals are given by: 
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Figure 1: Contours of integration to Eq. (1431) 



1 f°° e~ * 
I n( E,T) = -j_J x -^h(T, x) . (43) 

In these cases n = 9. These integrals can be evaluated choosing the contours 
given by Fig. (CQ) and using the residue theorem done in Refs. (TJ [18] . 
However, note that we have odd order poles. In fact, these integrals will 
be associated with a nine order residue of h(x) at % = 0. So, performing 
such integral the correction to the response function for the Unruh-DeWitt 
detector is: 



Wi(E,T) = ^^ xF(E,a) (44) 



where 



(-298262a 6 - 29302a 4 £ 2 + 12817a 2 £ 4 + 1842E 6 ) 
^ ,a '~ 12607ra ^ ' 

and now cr 2 (r) w cr 2 sinh(ar) w <7 2 ar, when ar 1. 

We note that the correction to the response function shows a Fermi- 
Dirac factor instead of a Bose-Einstein factor found in the free case. If 
a temperature could be defined here, it would be given by T = 2n/a. 
However there is some remarks to discussed here. Remember that in the 
Minkowski reference frame the noise function /x(x) depends only on the 
spatial coordinates, while in the Rindler system the noise function depends 
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on the proper time r, in other words the noise is a function of £ and 77, 
and now we have a non-equilibrium behavior of the field measured by the 
detector. In this case we can define a local temperature, see Refs. [21] and 
[22] . which depends on the proper time by the function 77). It was showed 
that the temperature which appears in the response function correction is an 
averaged temperature over the fluctuations and the global temperature is 
related with the thermal equilibrium. 

The last result shows that the lack of stationarity in the accelerated 
frame does not necessarily implies that a thermal-like response function with 
local temperature is forbidden. With such above considerations a possible 
interpretation of our result is that light-cone fluctuations in this model lead 
to the emergence of a local temperature as seen by accelerated observers. 

So far we have evaluated the horizon fluctuating contribution to the 
response function for an Unruh-DeWitt detector coupled to a massless 
fermionic field. Instead using a detector like that one considered in Eq. 
(|30p . it is also possible to use another kind of detectors. It is known that 
the interaction Hamiltonian for the scalar field coupled with the Unruh- 
DeWitt detector is a monopole- moment [12] . In other words, it is a two level 
system operator coupled with a scalar field, hence it is linear in the field. 
It gives a response function which depends on the two point instead of the 
four point Green function found to the fermionic field case. Thus, instead of 
Hamiltonian given by Eq. f!3Up . we can find an analogue kind of monopole- 
moment detector similar to the scalar case Hamiltonian. Despite the fact 
that this detector is less realistic, one can consider it equipped with a spinor 
and coupled linearly to the field via the Hamiltonian: 

H int = 6(r)*(z(r)) + ¥(x(t))0(t), (46) 

where x M (r) is again the world line of the two level system with respect to 
the proper time r, and 0(r) is the fermionic monopole moment operator. 
To keep the Hamiltonian stationary to the observer's proper reference frame 
it is also necessary to Fermi- Walker transport the monopole operator and 
the Dirac field, in same way that was done above using Eqs. ( I23i) -( l25j) . One 
can also formulate a perturbation theory in the interaction picture where it 
is assumed that = S T Q and = QS T rather than and obey the 
time-evolution equations. 

Following these steps one can find the new response function: 
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7° 



f dr' [ T dr"e- iE ^- T "^ (47) 
Jo Jo 

x(O|^(x(r'))^(x(r'0)|O)" 



This result is similar to that one found in Ref. [23] . 

The response function is proportional to the two point Wightman function 
instead of the four point function found previously. The calculation is done 
in a finite time interval as was done in Eqs. (14"T| |4"2|) . and up to first order in 
r we have that the correction due to horizon fluctuations is 

^ (-8059^ -55^ + 3684^) 
v ' ; 60vra(e 2 ^/ Q - 1) V ' 

the cr 2 (r) is the same as in Eq. (1441) . This integral was evaluated using Eq. 
(143]) with n = 6 and we also considered the same contour given in Fig. ([T|). 

In the case without fluctuations the two point function is given by Eq. 
(1281) . and the response is 

W a{ E,r) ^Kl {a f\ (49) 
e « +1 

Comparing Eqs. (148|) and (149|) we note that the event horizon fluctuation 
modifies the Fermi-Dirac response found in the free case. It appears a 
Bose-Einstein factor. As we discussed, it is also possible to define a local 
temperature, T = 27i/a, in this non-equilibrium situation and, in a similar 
way, the effective noise strength is proportional to the proper time. 



5 Conclusions 

Recently it was proposed in condensed matter physics an analog model for 
quantum gravity effects [I] and in Ref. [2] the authors extended this theory 
to a real massive scalar field. In Ref. pQ it was considered a massless scalar 
field near a four- dimensional Schwarzschild black hole with event horizon 
fluctuations. They used a perturbation theory similar to that one previously 
considered in fluctuating disordered media. 

In this paper we have used these ideas to study the behavior of a massless 
fermionic field near a fluctuating black hole horizon using the techniques 
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described before. First, we developed a perturbation theory associated to a 
massless fermionic field in the presence of random noise. After performing the 
calculations we have applied these results to investigate the thermal radiation 
near a fluctuating horizon. We showed that horizon fluctuation implies that 
the Unruh-DeWitt detector which has a Bose- Einstein factor in the response 
function associated with the free case has now a correction which has Fermi- 
Dirac factor. We would like to stress that we are assuming that it is possible 
to define a local temperature in this non-equilibrium situation. In a less 
realistic detector, coupled linearly with a Dirac field, the correction has also 
an apparent change in the distribution measured by in the response function 
with unchanged temperature. In this situation the leading term (case without 
fluctuations) has a Fermi-Dirac factor and in the correction appears a Bose- 
Einstein factor. 

Since the de Sitter horizon has also a thermal spectrum measured by 
the Unruh-DeWitt detector [21] for both, bosons and fermions, a natural 
extension of our work is to study the effects of a fluctuating de Sitter horizon 
to such kind of coupling. It is also natural to enquire the effects of introducing 
a time dependent randomness [25] in place of the Eqs. (J2j) and (JSJ. Such 
random coefficients allows the possibility of particle production from the 
beginning [TJ. These subjects are under investigation by the authors. 

Appendix A One-loop correction to positive 

frequency Wightman function 

Let us begin from Eq. f JT7|) which gives the correction to the positive 
frequency Wightman function. The self energy is given in Eq. (fT8j) . 
By dimensional regularization, only the first term do contribute which 
corresponds to 2itiuj 2 ^ Q . Replacing this result in Eq. ( 1T7|) this expression can 
be separated in three different parts ■y^'j 5 k fJi ks = SijU 2 — 2ukj^^° — Yl^hkj 
writing Eq. ( Tl6|) as 



gW(t,x) =2<K 2 ia 2 i 



o 




where I n are integrals defined by: 




(51) 
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To calculate the positive Wightman function, we must evaluate Eq. (1511) 
over the contour of integration for the variable u as is defined in Ref. [2E| . 
In this case only the pole lu = |k| is chosen and the residues theorem is used. 
It yields to 

4 = ^ n - 3 [(n - 1) - i\k\At]e- Wt . (52) 

The expression above is simplified and can be written in terms of the 
momenta as 



x) = ^ i^L HAf|k| + 5] 

Ikl 2 

+ 20^^ -j- [-iAt\k\ + 4] 

-lY^kikj^ [-iAt\k\ + 3] | e -*MAt+ik.Ax (53) 

The last two integrals can be simplified using the fact that, 
/ rf 3 k^77(|k|)e- ifc(x - x,) = fdij d 3 kf(\k\)e- ik{x - x '\ So our integrals are 
from the type 

B [ n) = J ^|kre- ife (*-*'), (54) 
the expression is simplified by 

S«(t,x) =27r 2 za 2 7° {^Bf ] + ±AtB[ A) - 2 1 ^°d,B{ 2) 

+ l - 1 i 1 °Atd j Bf ) + jj^didjB^ + ^yd^A . (55) 

After calculating the angular parts from the above integrals it is useful 
to define a function g n (x,x') such that the function B^ and its derivatives 
can be expressed in terms of g n (x, x') by the following expression 

poo 

g n {x,x')= dkk n (e- ik( - At ~ A ^ _ e - ife ( A * +A W)). (56) 
Jo 
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In the last expression k is a real variable not a four-vector, thus, it can be 
put out from the integral through derivatives with respect to At writing the 
last equation as: 



9n(x,x') 



I an poo 

This expression is related to the Wightman positive function for the scalar 
field, in such a way that the g n (x,x') function may be written by (see Ref. 
[26]) 

1 d n f -2i|Ax| \ 
9n{x, x ) - TZ^MF ^(At-ze) 2 -|Ax| 2 ) ' (58) 

Thus JBj 's integrals can be simplified by the expression 

Finally, substituting Eqs. and (J55|) into Eq.(l55j) we get the first 

contribution to the Wightman function which is given by Eq. (jl9p . 



Appendix B One-loop correction to the 

Rindler noise 

With the Rindler coordinates defined in Eqs. ff22l we are able to define the 
Rindler noise for the Dirac field by: #1/2 (t) = S t (S^(x(t),x(t')))S^ t /. So 
using this definition together to Eqs. ()19p . fl20|) and (12 ip we find that in 
accelerated reference frame 



2D G 



7°^ - 



a 



a 



2D 



a 



sinh 5 — 
2 

sinh 5 — 
2 



"2 at 

-smh^F(z) + ieG(z) 
a 2 

2 at; 

— sinh — if (z) + ieJ(z) 
a 2 



(60) 



where we have defined z 
defined by 



The functions F(z),G(z), H(z) and J(z) are 
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F(z) = cosh z (—3684 cosh 5 z + 7539 cosh 3 z sinh 2 z 

-8328 cosh z sinh 4 z - 1005 cosh" 1 z sinh 6 z) , (61) 

G(z) = cosh z (18504 cosh 4 z - 24732 cosh 2 z sinh 2 z 

+9354 sinh 4 z) , (62) 

i?(z) = cosh z (—1660 cosh 4 z sinh z — 2024 cosh 2 z sinh 3 z 

- 156 sinh 5 z) , (63) 

J(z) = cosh z (6620 cosh 3 z sinh z + 5056 cosh z sinh 3 z 

+ 12 sinh 5 z cosh" 1 z) , (64) 

expanding the denominator in Eq. (1601 and after some calculations we find a 
simplified expression for the correction of the Rindler's noise given in Eq. (l29|) . 
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